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Abstract 

This paper proposes lower bounds on a quantity called ZA-norm joint spectral 
radius, or in short, p-radius, of a finite set of matrices. Despite its wide 
range of applications to, for example, stability analysis of switched linear 
systems and the equilibrium analysis of switched linear economical models, 
algorithms for computing the p-radius are only available in a very limited 
number of particular cases. The proposed lower bounds are given as the 
spectral radius of an average of the given matrices weighted via Kronecker 
products and do not place any requirements on the set of matrices. We show 
that the proposed lower bounds theoretically extend and also can practically 
improve the existing lower bounds. A Markovian extension of the proposed 
lower bounds is also presented. 

Keywords: p-radius, switched linear systems, mean stability, Markov 
processes 


1. Introduction 

The L p -norm joint spectral radius [1, 2] (often called p-radius ) of an in¬ 
dexed family of n x n real matrices M. = {Ai,..., A^v} is defined by, for 
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where ||-|| denotes the maximum singular value of a matrix. Introduced by 

Jia [3] and Wang [4] independently for p — 1, the p-radius now plays an im¬ 
portant role in various fields of applied mathematics. A classical application 
of the p-radius is in the characterization of the regularity of wavelet functions 
in L p spaces [4, 1], The mean stability of a class of switched linear systems is 
determined by the value of the p-radius [5, 6]. The so-called indeterminacy of 
a switching linear economic model can be checked through the p-radius [7]. 

However, being defined as a limit of a sequence on matrix products, the 
p-radius is known to be difficult to compute. There is no formula available 
for its computation except under the conditions that p is an even integer 
or that the given set of matrices leaves a common proper cone invariant [2, 
5]. Although for the latter case there exists a converging approximation 
method [8] that does not require p to be an integer, in a general case, even 
approximating the p-radius is an NP-hard problem [8]. As for bounds on 
the p-radius, the sequence defining the p-radius is decreasing and therefore 
gives upper bounds, though their computation requires exponentially growing 
costs. Finally, the lower bounds on the p-radius in the literature [9, 7] are 
often not very accurate. 

In this paper we propose novel lower bounds on the p-radius for integer 
values of p with no assumptions on the given set of matrices on the contrary 
to [2, 8]. The lower bounds are given as the spectral radius of a weighted 
average of the given matrices and the weights are realized by Kronecker 
product of matrices. We will show that, with appropriately chosen weighting 
matrices, the proposed bounds extend and also can improve the lower bounds 
in the literature [9, 7]. The obtained results are furthermore generalized to 
the Markovian setting. This in particular enables us to use the p-radius 
to study the stability of so-called Markov jump linear systems [10], which 
are switched linear systems whose parameter changes by following a time- 
homogeneous Markov chain. The generalization is based on a stochastic 
counterpart of the so-called O-lift of matrices [11]. 

The organization of this paper is as follows. In Section 2, we derive a 
novel lower bound for the p-radius. In Section 3, we provide a Markovian 
extension. The notations used in this paper are standard. The identity 
matrix is denoted by I. The spectral radius of a square matrix is denoted 


2 



by p(-). The Kronecker product (see, e.g., [12]) of matrices A and B is 
denoted by A <8> B. Ai denotes an indexed family { Ai ,..., ofnxn real 
matrices. Finally, for ease of reference, we list some important properties of 
the p-radius in the following proposition: 

Proposition 1. Let p > 1 be a positive integer. 

1. The sequence {hk{Ai)}^ =l defined by 



is decreasing [8]. 

2. If either p is even or Ai leaves a proper cone invariant, then p p (Ai) = 



where A® p denotes the Kronecker product of p 


copies of A (see [2]). 


3. It holds that p p (Ai ,..., A N ) = (pi(Af p ,..., A® p ))^ p (see [8]). 


Remark 2. Throughout the paper, we often omit the dependence of quantities 
on the underlying family Ai of matrices when it is clear from the context. 
Precisely speaking, for a function / defined on (M nxn ) Ar , we may simply 
write f(Ai) as /. Also, abusing notation, we sometimes write f(Ai) as 
f(A 1 ,..., A n ) when Ai = {A x ,..., A N }. 

2. Novel lower bounds 

In this section, we present novel lower bounds on the p-radius. We also 
prove that these new bounds outperform existing lower bounds. Notice that, 
according to the third claim of Proposition 1, any result on the 1-radius 
(p = 1) is directly applicable to the general p-radius; thus, we shall focus on 
the particular case p = 1 for the rest of the paper. In order to state our main 
results, we need to recall the definition of the joint spectral radius [13]: 


Poo (Ad) = limsup max 





Our first result is stated in the next theorem: 
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Theorem 3. For W = { W\ ,..., Wjv} C M mxm , let 


N 


X w {M) = p[ A 


1=1 


IfPoo(W) = 1, then \w(M) < pi(-Ad). 
Proof. Using Gelfand’s formula, we obtain 


Aw(.M) = lim 


/c—>-oo 
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Using general identities [12] about Kronecker products: 

(A <8> C)(B <g> D) = (AB) ® (CD), 

P®£|| = P|| PH, 

we can evaluate the norm in (1) as 
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Taking the power of 1/A; of the last expression, and substituting in (1), we 
obtain 

Aw (-Ad) < poo (W)MM) = (4) 

as desired. □ 

Remark 4. From (4), we can see that the equality constraint Poo(bV) = 1 in 
Theorem 3 can be relaxed to the inequality constraint Poo(W) < 1. Although 
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the largest lower bound is clearly attained when the joint spectral radius of 
W is equal to one, it is convenient in practice to relax the equality constraint 
because checking that Poo(kV) = 1 is NP-hard [13]. In contrast, checking the 
inequality Poo(W) < 1 can be efficiently done (using, for example, the JSR 
software toolbox [14]). 

In what follows, we compare the proposed lower bounds with the most 
relevant existing bounds found in the literature. The first lower bound is the 
one implicitly obtained by Zhou [9]; if we let 

“ - AT -> 

N Poo 

then pi > ii- This inequality can be derived from the following bound 
on the joint spectral radius: Poo > {p p+q /p p ) v ^ q p v+q whenever p,q > 1 [9, 
p. 48]. Letting p = q = 1 gives p\ > pi/ Poo- Then, applying assertion 2 of 
Proposition 1 to this inequality proves p\ > £%. The second lower bound was 
introduced in [7]: for wi ,..., Wjy £ [—1,1] it holds that 


Pi > P 



In the following theorem, we show that the lower bound herein proposed 
extends the ones mentioned above: 


Theorem 5. For each m > 1 define 

£ m (M) = sup {Aw : W £ (r xm ) JV , 

The following statements are true: 

1. If m> n, then I m > Iz- 

2. Define 


= max p 
1Ul,...,W);ve[-l,l] 

Then < £ m for every m > 1. 



Poo(W) = 



!}• 


(5) 


( 6 ) 
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Proof. We begin by showing that 

[m<m'} => [£ m <£ m '\. (7) 

Let us take an arbitrary W C (j^ mxmS j N such that Poo(VV) = 1. For each 
i we define the block diagonal matrix Vj = diag(IFi, G xm and 

let V = {Vi,..., Vn}. Then, it holds that Ayy = Ay- Since Poo(V) = 1, this 
implies that Ayy < £ m i. Finally, taking the supremum with respect to W in 
the last inequality proves (7) 

Let us prove the theorem. It is clear that the normalized family of 
matrices M.q = {Ali/poo,..., AIat/Poo} C M nxn has a joint spectral radius 
equal to one. Therefore, from Theorem 3 it follows that £ n > A m 0 — 
p(lV _1 Y^i =i A <8) Af)/poo = li- This inequality and (7) prove the first claim 
in the theorem. To prove the second claim, we take an arbitrary family 
{uq, ..., tDjv} C [—1,1]. Since this family has a joint spectral radius less 
than or equal to one, we obtain p{N l Ylf=i w iAi) = A{ Wli ... iWjv } < i\ by 
Theorem 3 and Remark 4. Therefore, < £\ and hence (7) proves the 

second claim. We finally remark that the maximum in (6) exists because the 
set [—1,1]^ is compact and the function p(-) is continuous. □ 


In what follows, we illustrate our results with some examples. 


Example 6. Let IV > 1 be arbitrary. Consider the matrix family A4 = 
{Ai ,..., ALat+i}, where A { = NI and 


A 2 


An+i — R 


0 -1 

1 0 


( 8 ) 


Since \\R k \\ = 1 for every k > 0, we see that pi equals the 1-radius of the 
family of scalars {IV, 1,..., 1}, where element 1 has multiplicity N. The 
1-radius of this family equals 


P 


(N+ 1 + --- + 1A 

V N+l J 


2 N 
N + l 


by assertion 2 of Proposition 1. Therefore p i(A4) = ( 2N)/(N + 1). This 
value is attained by the proposed lower bound £2 because 

A { I.R. . K) = P f Wy (I ® (NI) + NR ® R)\ 

N 

= j^p(I®I + R®R) 

2N 

~ N + l' 
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On the other hand, a straightforward computation shows that = 1. 

Example 7. Let M. = {/, R, R 2 , J? 3 }, where i? is given in ( 8 ). Clearly, 
Pi = 1 because ||i? fc || = 1 for every k > 0. The proposed lower bound attains 
this exact value of the 1-radius as £ 2 = 1 because A/ R R 2 R 3 = 1 . On the 
other hand, we can show that £[~i,i] < a/2/2 as follows. Let Wi G [—1,1] 
(i = 1,..., 4) be arbitrary. Then we have 

_ Wil + w 2 R + iv 3 R 2 + uqi ? 3 _ 1 u -v 
4 4 v u ’ 

where u = w± — w 3 and v = w 2 — uq. Since \u\ < 2 and |n| < 2, we can see 
p(S ) < a/2/2 and therefore /-i,i] < \/ 2 / 2 . 

Although can largely improve other lower bounds in the literature, it 
is not easy to compute. The first reason is the non-convexity of the func¬ 
tion p (see, e.g., [15]). The other reason is that the set {W G (K mxm ) Ar ; 
Poo(VV) = 1} does not admit an appropriate parametrization due to the NP- 
hardness [13] of computing p^. For these reasons, we here propose using the 
set of matrix weights from the following set: 

£> m = {{DiLi}^ = i : Di is diagonal, ||ZA|| < 1 , and L,: is orthogonal} . 

Since the matrices in this set have a joint spectral radius less than or equal 
to one, the best lower bound = max/Aw : W G O m } achieved by using 
such matrix weights provides a lower bound on the 1-radius by Theorem 3 
and Remark 4. Moreover, since an orthogonal matrix L admits [16] the 
parametrization L = D(I — S)(I + S ) _1 where S' is a skew-symmetric matrix 
and D is a diagonal matrix whose diagonals are either +1 or —1, we can max¬ 
imize £o(m) using, for example, a stochastic gradient descent algorithm [17]. 
Notice that the result of this algorithm will be a local maximum. 

In the following example, we illustrate the effectiveness of the weights 
from the set O m by studying the stability of a switched linear system. 

Example 8 . Consider the switched linear system 

X(k + 1) = A a{k+1) X(k), X( 0 ) = I (9) 

where {a(/c)}^L 1 are random variables independently and uniformly dis¬ 
tributed on {1,...,1V}. The system is said to be pth mean stable [18] if 
there exist C > 0 and 7 G [0,1) such that 

£[IW*0in < CY* ( 10 ) 
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for every k > 0. It is well known [13] that 

[(9) is pth mean stable] [p p < 1]. (11) 

In this example, we let IV = 2 and randomly choose two matrices 


'-0.87 

-0.77' 

, ^2 — 

'0.14 

0.40 ' 

1.17 

-1.09 

0.89 

-0.73 


Using a MATLAB implementation [19] of the stochastic gradient descent 
algorithm, we perform the maximization over the set O m and find 


'-0.71 

-0.70' 

, Wo = 

'0.85 

-0.53' 

0.70 

-0.71 

5 v v Z 

0.53 

-0.85 


Since \{w 1: w 2 } — 1-07 we obtain l 2 > 1.07. Tims, Theorem 3 implies that 
p i > 1 and hence the the system in (9) is unstable, according to (11). We 
cannot prove that the system in (9) is unstable using the other lower bounds 
in the literature, since lz = 0.93 and = 0.73. We remark that the joint 
spectral radius appearing in ti is evaluated with the JSR Toolbox [14]. Also, 
the maximum in U_i,i] has been evaluated with extensive simulations on the 
weights. 

We propose a further extension of our bounds based on product families 
of a set of matrices. Let us define A4 q to be the family of matrices consisting 
of all the N q products of the matrices from M. having length q. Then, our 
extension can be stated as follows: 

Theorem 9. Let m and q be positive integers. Define £m (M) = £ m (M q ) 1/q . 
Then, £m < pi. Furthermore, if q is a divisor of another positive integer q', 
then £m < £m\ 

Proof. Let us first recall the following identities [13]: 

Pr(M q ) = Pp (M)\ (12) 

Poo(M q ) = Poo (M) q . (13) 

Then, using Theorem 3 and Equation (12), we can prove the first claim in 
the theorem as pfiM) = pi(M. q ) 1 ^ q > £ m {M q ) l ^ q = £m{M). To prove the 
second claim in the theorem, we let q = 1 and q' — 2 for simplicity. The proof 










for general q and q' is similar to this particular case and hence is omitted. Let 
W = {Wi ,..., W N } C M mxm be arbitrary and assume Poo(W) = 1. Since 
p(M) 2 = p(M 2 ) for a square matrix M, we can show that 


\w(M ) 2 



A MM 2 ). 


(14) 


Also, since Poo(VV 2 ) = Poo(hV) 2 = 1 by (13), it follows that A ^(M 2 ) < 
£(M 2 ) = ^ 2 \M) 2 . This inequality and (14) yield Aw < f^ 2 A Taking the 
supremum with respect to W in the left-hand side of this inequality proves 
iin < £rn , as desired. 


□ 


We close this section by giving a remark on complex weights. 

Remark 10. In principle, one could obtain better lower bounds using complex 
weights in £ m instead of the real weights W±, ..., Wn- However, we can show 
that it does not lead to an essential improvement. Precisely speaking, we 
here prove the following claim: if W = {Wi,...,Wjv} C C mxm satisfies 
Poo(W) = 1, then 

Aw A ^2 m- (15) 

To prove this claim, for W G C mxm we let 


T w — 


ReW 

Irn W 


- Im W 
ReW 


e ]R( 2m ) x ( 2m ) 


where Re W and Im W denotes the real and imaginary parts of W, re¬ 
spectively. The multiplicative property Tww = T\yTw' and the identity 
\\Tw || = \\W\\ yield 

Poo (T Wi ,...,T W n ) = 1. (16) 
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Also, since p(W) = p(T w ), we can show that 


Ayv(.M) 


P ( T n - 


N - 1 E i=i Wi®A, 
N 


P 


P 


Ef=i(ReWi)® A 


N LE^i(ImWi)® A 


TV 


— ( 8 ) 4 


2 = 1 


AjT^ ,...,Tw n } (-A^) • 


EL-Omw;)® 4l\ 

Eli(ReWi)® A Jy 


This equation and (16) prove Inequality (15). 


3. The Markovian case 

In this section, we extend the results presented in the last section to the 
Markovian case. Let o = {cr(/c )}^ =1 be a time-homogeneous Markov chain 
with a state space {1,..., N} and a transition probability matrix Q e R JVxAr . 
We define the Markovian version of the p-radius as follows. 

Definition 11. Let X(-; p) denote the trajectory of the Markov jump linear 
system 

E : X(k + 1) = A a{k+1) X(k), X(0) = /, cr(l) ~ p, 

where p is an arbitrary probability distribution on {1,..., N}. The L p -norm 
Markovian joint spectral radius ( Markovian p-radius for short) of the pair 
(M., Q) is defined by 

p p (M,Q) = sup limsup [E[\\X(k] /i)|| p ] 1/(pfc) ) . 

fi k —^oo 

Remark 12. The logarithm of the Markovian p-radius corresponds to a quan¬ 
tity called Lyapunov exponent of the pth mean [20, p. 307]. 

We will later see that this definition coincides with the originally defined 
p-radius in the case where the matrices are independent and uniformly 
distributed at every time step. Moreover, as can be naturally expected, the 
Markovian p-radius has a close connection with the mean stability of E, which 
is defined as follows. We say that E is pth mean stable if there exist C > 0 
and 7 e [0,1) such that (10) holds for all k and p. The next proposition is 
an immediate consequence from the definition of the Markovian p-radius. 
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Proposition 13. £ is pth mean stable if and only if p p (A4,Q) < 1. 

Proof. The necessity is obvious. Let us prove the sufficiency. Assume that 
p p (M., Q ) < 1. Then there exists a 7 G [0,1) such that, for every i G {1,..., N}, 
we have limsup fc ^. 00 (if[||X(/c; 5j)|| p ] 1/, ( pfc ^) < 7 , where <5 i is the probability 
distribution on {l,...,jV} such that 5j({i}) = 1. Then, for each i there 
exists a positive integer AT* such that, if k > Ki, then E[\\X[k\ 5«)|| p ] < 7 pfc . 
Therefore, for an arbitrary /q if k > max (Ah,..., Kn) then E[\\X(k] /i)|| p ] = 
YliLi PiE[\\X(k] 5j)|| p ] < YliLi lMl pk = l pk - This implies that £ is pth mean 
stable. □ 


The main result in this section is stated in the following theorem, which 
can be used to compute upper and lower bounds on the Markovian p-radius. 

Theorem 14. 1. The sequence {h k (M., <5)}fcT 1 defined by 

( N \ i/(pfc) 

^ ^ Qh,i2 ' ' ' Qik-i,ik II^4 ' ' ' l| P j 


is decreasing. Moreover p p (Xt, Q ) = Hindoo h k (A4 , Q ). 
2. For VV = {Wij} \<i.j<N C R mxm define 


— 


Aibbn ® A 1 

qwWiN <S> Ai 


QniW N i ® A n 
QnnWnn <S> A n 


If PoJW) = l, then pi(M, Q) > p{A\^). 

Let us observe some consequences of Theorem 14 before proving it. First, 
consider the special case when a is a sequence of independent and uni¬ 
formly distributed random variables on {1,..., N}. The corresponding tran¬ 
sition probability matrix Q is the N x N matrix whose entries are all 1/N. 
In this case, we have h k (M,Q ) ={Y^ 1 ....,i k =i{ l / N ) k ~ l \\A ik - ■ ■ A^) 1 /^ = 
N l t k Ph k (M). This equation implies that h k (Xi,Q), the average of the norm 
of /c-product A ik ■ ■ ■ A ix , coincides with the other average h k (M) under the 
identification of £ as a switched linear system having independent and iden¬ 
tically distributed system parameters, except the factor A rl//fcp . This factor, 
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roughly speaking, arises because h k (M.,Q ) does not take the initial prob¬ 
ability distribution of the switching signal a into account. Then, by this 
equation and the first claim of Theorem 14, taking the limit as k —> oo 
shows that p p (M.,Q) = p p (Ai). Hence Definition 11 indeed recovers the 
original p-radius. 

Next, as a corollary of Theorem 14, we can recover a lower bound of the 
Markovian p-radius implicitly presented in [7]. 

Corollary 15 ([7]). Define 

tm(M, Q ) = sup |p(Ay) : W G (r xm ) Af2 , Poo (W) = l} . 


Then, £i(M, Q) < pi(A4, Q ). 

Before we prove Theorem 14, we illustrate its use with an example. 
Example 16. Consider the Markov jump linear system £ with 


Tli 


" 0.77 

0.80' 


'-0.77 

0.83 ' 

-0.60 

0.87 

, a 2 — 

-0.70 

-0.70 


Q 


0.70 0.30 
0.43 0.57 


A brute force search indicates that £i(M, Q ) = 0.844 with the weights W n = 
1, W 12 = 1, W 2 1 = —1, and W 22 = 0.932. On the other hand, using a 
MATLAB implementation [19] of the stochastic gradient descent algorithm 
used in Example 8, we locally maximize 0(A4, Q) over the set D 2 of matrix 
weights to find the set W consisting of the following matrices: 


'-0.412 

-0.911' 

, w 12 = 

'0.839 

-0.544' 

0.911 

-0.412 

0.544 

-0.839 

'-0.204 

-0.979' 

, w 22 = 

'0.937 

-0.349' 

0.979 

-0.204 

0.349 

-0.937 


Since p(Ayj) = 1.067, we conclude pi(A4,Q) > 1.067. This proves that the 
corresponding Markov jump linear system is unstable. 

The rest of this section is devoted to the proof of Theorem 14. In our 
proof, we employ a reduction of the Markovian p-radius to the original p- 
radius, for which we can apply the results obtained in Section 2. 
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Proposition 17. Let p be a positive integer. For each 1 < i < N, let e* 
denote the ith vector in the canonical basis of~R N . Define the set of matrices 
Ai p {Bij 


B$> =N 2 'rqt/’’A i ®(e j el). 


Then, 

Pp(M, Q) = Pp(Mp). (17) 

Proof. We first claim that, for the proof of (17), it is sufficient to show that 

\p P (M,Q)< 1] \Pp(M p )< 1], (18) 

due to the following reason. Suppose that (17) does not hold while (18) is 
true. Then, we have either p p (Ai,Q) > p p (A4 p ) or p p (A4,Q) < p p (A4 p ). If 
the former inequality holds, then one can find a c > 0 such that the matrix 
family cAi = {cA i } I f =1 satisfies p p (cAi,Q ) > 1 > p p (cA4 p ) = p p {cAi p ). 
However, this cannot be true by (18). In a similar way, we can also show 
that the latter inequality, p p (M,Q) < p p (M. p ), cannot hold. Therefore, (17) 
must be true. 

Consequently, it suffices to prove (18). In order to prove the claim in (18), 
we introduce an alternative switched linear system with independent and 
identically distributed jumping parameters. Let {6 ) (/c)}^ =1 and {4>(k)}kLi be 
independent random variables uniformly distributed on {1,..., IV}. Define 
the switched linear system by 

S p : X(k+ 1) = X(0 ) = L 

From (11), we see that p p (Ai p ) < 1 if and only if is pth mean stable. 
Moreover, from Proposition 13, we know that p p {Ai, Q) < 1 if and only if £ 
is pth mean stable. Therefore, to prove (18), we need to show that £ is pt\i 
mean stable if and only if £ p is pth mean stable. 

To prove the equivalence of stability, let us first compute £'[||W(A;)|| P ]. By 
the deffiiition of £ p and Equation (2), we can compute X(k) as 

k 

m =n (^C®^ ® (^w4 (i) )) 

i =1 

= N 2k / p (q d { i),</>(i) • •' Qe(k),ct>(k)) ^ (^e(fc) • ■ • ^4e(i)) ® Jk , 
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where the symbol n^Li denotes the left product of matrices and also = 
e </>(fc) e e(fc) ''' e </>(i) e e(i)- Since the vectors ei,..., e# are orthonormal, ||J fc || = 1 

if 

= 9(i + 1), i = 1,..., k - 1, (19) 

and otherwise ||Jfc|| = 0. Therefore, if we denote by x the characteristic 
function of the event (19), using (3) we obtain 

||X(fc)|| p = xN 2k qe{i)m ■ ■ • qe{k-\),e{k)qo{k),<t>{k)\\Ao(k) ■ ■ • -4<?(i)|| p - 

Since the event in (19) occurs with probability l/iV fc-1 , we obtain 

E[\\X{kW] 

N 2k 

= J^ 1 E [^(1),0(2) • • ■ qe(k-i),e(k)qe(k)^{k)\\Ae(k) ■ ■ ■ A 0 (i)|| p ] 

N 1 

^ ^ 1 jyk+l J 11-^-ifc ' ' ’ l| P 

ii,.. .,4,7=1 (20) 

N 

— 'y ] 9u,t2 ■ ' ' ?4-i,4 11^4 ' ' ' -^ii II 

^1 VJ^fe = l 

N 

= ^E[iiA'(*i*,)r], 

H —1 

where we used X!)=i ? 4 ,j = 1 to show the third equality. 

Now, assume that £ is pth mean stable. Then, by (20) we have i7[|| A - (k) || p ] 
< Cx pk = CN^ pk for some C > 0 and 7 G [0,1), and hence £ p is pth 

mean stable. On the other hand, assume that £ p is pth mean stable. Then, 
there exist C > 0 and 7 G [0,1) such that £'[||X(/c)|| p ] < C' 7 pfc . Therefore, 
(20) shows that E[\\X(k; < s >i 1 )|| p ] < C~/ pk for every i\ G {1,..., N}. Then, in 
the same way as the proof of Proposition 13, we can conclude the pth mean 
stability of £. This completes the proof of the theorem. □ 

Remark 18. From Proposition 17, we can regard the matrix B ( - p> as an L p - 
averaged version of the O-lift introduced in [ 11 ], which is used to generalize 
the so-called Berger-Wang formula to a Markovian version of the joint spec¬ 
tral radius. Also, we remark that considering the auxiliary switched linear 
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system E p with an extended state space is similar to considering an ex¬ 
tended state-variable consisting of the original state variable and the under¬ 
lying Markov chain, which is frequently employed for studying (semi-)Markov 
jump linear systems [10] ([21]). 

Finally, let us prove Theorem 14. 

Proof of Theorem If. From Equation (20), we see that 

h k (M p ) = E[\\X(k)r\ l /W = h k (M,Q). 

Therefore, by Proposition 1, the sequence {h k (A4 } Q)}^L 1 is decreasing. Fur¬ 
thermore, it converges to lim^oo h k (M. p ) = p p (Al p ) = p p (Ai,Q) by (17). 
Thus, the first claim is proved. Let us then prove the second claim. Assume 
that VV = {W t j}\<i.j<N C M mxm satisfies Poo(VV) = 1. By (17), it is enough 
to show that pi(Ali) > p{Ayf). From Theorem 3 we obtain 


Pi(Mi) > p 


= P 



N 

E w » ® df 


N 


Y Qij W ij ® A i 


\i,j =1 



( 21 ) 


Here we recall that there exists [12] an invertible matrix T satisfying C®D = 
T~ 1 (D <g) C)T for all C G ~R NxN and D G x ( nm ). Consequently, there 

exists an invertible matrix T such that W t j <g) Ai <g> ( ejej ) = T^^ejeJ) <8) 
Wij <g) Af)T. Therefore, the matrix appearing in the last term of (21) is 
similar to Yij=i Qij( e j e J) ® ® A, which in fact equals Aw- Therefore 

Pi (All) > p(Mvv), as desired. □ 


4. Conclusion 

This paper proposed novel lower bounds on the p-radius of a finite set 
of matrices. The obtained lower bound is given by the spectral radius of an 
average of the given matrices weighted via Kronecker products. We showed 
that the proposed lower bounds theoretically extend and also practically 
improve the existing lower bounds. We have also shown the extension of the 
p-radius and its lower bounds to the Markovian case. 
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